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This is the final report describing the resampling parfcrscd by 

IBM for the NASA Goddard Space Flight Center under modifications 6 
through 9 of contract NASS-21865. The objective of this effort was to 
compare the nearest neighbor and cubic convolution resampling algorithms 
as applied to Landsat MSS data. The conq>arlson was made by evaluating 
the effects of the resampling algorithms on each of the following cate- 
gories : 

o "Point source" objects 

o Spatial frequency content of the data 

o Areas of uniform radiance 

o Multlspectral classification 

o Temporal registration 

o Human analysis of the images 

This report describes the work that was performed, presents copies of 
the results generated, analyzes those results, and presents conclusions 
drawn from them. 
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SUMMARY OF EFFORT 


Five groups of data seta were procesaed for this Investigation. Each 
was resampled with the nearest neighbor end cub'c convolution resampling 
algorltljDS , and a comparison of the resulting dat-i arrays was made by 
evaluating various photographic and statistical proaucts. For a given 
data set, the sampling lattice used for nearest neighbor resampling 
was Identical to that used for cubic convolution. Any differences in 
the resulting digital arrays are therefore due solely to the properties 
of the resampling algorithms. In preparing the photographic products, 
the original and resampled versions of each data array were processed 
Identically. Hence any differences In the resulting prints are also 
due solely to the resampling algorithms. 

All film products were generated from recordings made on the IBM Drum 
Scanner/Plotter using a 100 pm square spot. For data arrays whose 
pixels have a 1.4:1 aspect ratio, the ground scales of the film 
recordings are approxlsutely 1:560000 horizontally and 1:790000 
vertically. For arrays whose pixels have a 1:1 aspect ratio, the 
ground scale of the recordings Is approximately 1:500000 In both 
directions. Ten-power enlargements were made of all film recordings. 
These have correspondingly Increased ground scales (e.g., 1:79000). 

The data sets processed are listed In Figure 1. The generated products 
are listed In Figure 2 and described below. 

Data Set A. - Point Source Objects 

To determine the effects of the resampling algorithms on single-pixel, 
bright objects, a data set containing Images of the mirrors placed by 
U. E. Evans In the vicinity of Reston, Virginia in October 1974, was 
used. The Intensity values surroui.dlng each mirror In the original 
and resampled arrays were printed. Since the mirrors were within 
four or five pixels uf each ether, an 18 pixel by IS line array was 
sufficient to contain all th« -eievant data values. These array print- 
outs are presented In Flgurei. 6 through 10. Flgtre 3 contains a contact 
print of the .ata set, and Figures 5 threu^i 7 contain 10-power 
enlargements. The locations of the mlii.'r>^ are marked on the photo- 
graphic products and data printouts. 

Data Set B. •• Spatial Frequency Content: of the Data 

Three arrays were selected to contain features with high spatial frequency 
content, such as road networks and airports. Each array, in both its 
original and resampled forms, consists of 128 lines of 128 pixels per 
line. To determine the effect of resampling on spatial frequency in 
the horlzonte'i direction, each row (horizontal line) in each array was 
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trantforaed using the Fast Fourier Transfona (FFT) algorithm. For each 
line* the power In each frequency was computed from the transform, and, 
for a given array » the average power was computed as the component -wise 
average of the power of the 128 lines. To determine the comparable 
effect In the vertical direction, this process was repeated for each 
column (vertical line) the three arrays. Three standard operations 
were applied to each line of the data before It was transformed: 

(1) The mean was removed. 

(2) A window was applied (the window was a cosine bell applied 
to 1/10 of the data at each end) . 

(3) 128 pad values of zero were appended. 

The resulting power plots are shown in Figures 20-37. 

In addition to the spatial frequency analysis, the original and resampled 
arrays were plotted, contact printed (Figure 4), and enlarged (Figures 11 
19). The generated photographic products were evalued by an experienced 
photointerpreter as to overall quality, clarity, resolution, and infor- 
mation extraction potential. 

Data Set C. - Areas of Uniform Radiance 


To determine how the resampling algorithms affect areas of near uniform 
radiance, two such areas were selected: one near Lake of the Woods, 

Virginia, and the other In the Chesapeake Bay southeast of Point Look- 
out. The two areas were resampled with both algorithms. Shadeprints 
are presented in Figures 38 through 43. 

Data Set D. - Multlspectral Classification and Temporal Registration 

A systematic geometric correction was defined for a subimage of scene 
2183-16433 which contained Lade Intensive Test Site 92 in Hand 
County, South Dakota. Geometric correction parameters were then derived 
which would register the corresponding area of scene 2165-16435 with 
It. The Hand County site was chosen for the availability of detailed 
ground truth maps for a 5 x 6 nautical mile segment . The ground cover 
in this ragment, with the exception of about 2 % of the area, consisted 
of spring wheat, oats, com, and pasture/grasses, their exact pro- 
portions being known from the ground truth Inventory. 

The ground truth Information available for this site Is included as 
Appendix A of this report. It shows that the fields involved varied 
In size from 1 to 635 acres. Of the 275 fields, 60% were 40 acres or 
less In size, 27% were from 40 to 100 acres, and 13% were larger than 
100 acres. Approximately 55% of the fields had aspect ratios of 2:1 
or greater. In most cases, the field boundaries were regular and 
aligned along section lines. 
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All four bands of aach sublmaga wore raaamplad with the two resampling 
algorithms. Contact prints (Figures 44 and 45) and lOX enlargements 
(Figures 46 -69) for all four bands of each of the two temporally separ- 
ated scenes were generated and evaluated by an experienced photo interpreter. 
(Notice that Figures 50 - 57 and 62 - 69 are rotated to fit on the page, 
even though the figure titles are not.) 

The temporal registration that was done for this data set, and for data 
set e below, used relative control points that were correlateu between the 
reference image and the register image. Prior to the correlation, the proto- 
type features (the windows) and the search areas were resampled with a 
one-dimensional cubic convolution algorithm to remove the effect of sampling 
delay and intersweep earth rotation. Based on the locations of the window 
areas within the corresponding search areas, geometric correction parameters 
which would effect the temporal registration were derived. This one set of 
parameters was used to drive both the nearest neighbor and cubic convolution 
resampling algorithms. 

The digital arrays generated by the resamplings were classified using the 
multlspectral classification capabilities of the Earth Resources Interactive 
Processing System (ERIPS) at IBM - Houston. The accuracy of the proportion 
estimates obtained by the classification of each resampled version of the 
data was used as a basis for evaluating the resampling techniques. 

The classification model which was used in this study assumes that the 
data corresponding to each distinct ground cover constitute a normal 
population. This Is the model generally adopted In the classification 
of remotely sensed data. Implementation of the maximum likelihood 
classification scheme with this model consisted of the following steps: 

o Designation of training areas for each distinct ground cover from 
the ground truth map? 

o Computation of the statistics (mean vectors and covariance 
matrices) to characterize each type of ground cover (class) . 

o Classification of the data into the defined classes by a maximum 
likelihood rule. 

o Generation of classification maps and/or summaries of classlfl- 
results for training and test fields. 

There were mainly two considerations in the selection of training fields 
from the ground truth Inventory. First, It was desired to define training 
areas to be entirely in the interior of fields so that the boundary pixels, 
where the larger differences due to resampling show up, would not be Included. 
This eliminated narrow, small fields from being Included in training. The 
second consideration was to give a roughly proportional representation in 
terms of size and number of training fields to each of the four major ground 
cover types: spring wheat, oats, com, and pasture/grasses. Once the 

training fields were defined, mean and standard deviation statistics were 
computed for each class. Using the maximum likelihood decision rule, the 
test and training fields were classified for each image separately . Each 
pixel was assigned to one of the classes defined by the training fields. 
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Th« artal proportion estlmnt* of each grotmd cover category was computed 
directly aa the proportion of plxela asalgned to classes defined by the 
training fields of this category. As noted earlier, roughly 22 of the area 
had ground cover unlike any represented by the training fields. It was 
assumed that these pixels would be classified luto the defined categories 
proportionally by area. The proportion estimates obtained are stismarlzed 
In Figure 70. 

Data Set E. Temporal Registration 

Scene 1921-18022 was used as the reference Image and 1813-18063 was the 
register Image. Only band 5 of two subimage areas was processed. Contact 
prints (Figures 80 and 81) , lOX enlargements (Figures 82-93) , and shade- 
prints (Figures 94-lOS) were generated. The photographic products were 
examined by an experienced photointerpreter. 
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Section 2. EVALUATION AND CONCLUSIONS 


The output generated by the effort described above was evaluated* resulting 
In the following observations and conclusions: 

1. Cubic convolution may Introduce **overshoot” Into the resampled 
data. This effect Is clearly seen In the black borders sur^ 
rounding the mirrors In the vicinity of Res ton, VA. (Figures 7A 
and 7B) and the small white area at the left center of area 2, 
1921*18022 (Figure 87). Although not demonstrated In this In- 
vestigation, theoretical considerations Imply that the overshoot 
Is highly dependent on the sampling lattice. A slight shift of 
the input space relative to the sampling lattice may produce 
significantly altered resampled values at the same lattice 
points* In this sense, and In this admittedly extreme situation, 
cubic convolution is an Inaccurate algorithm for interpolating 

to inter-pixel values. 

2. Cubic convolution may cause a significant spreading of small 
features. This spreading was noticed more severely in the 
vertical (along-track) as opposed to the horizontal direction. 

In the original and nearest neighbor data, the southeastern 
mirror near Reston, VA appears to have affected only one hori- 
zontal line, but in the cubic convolution data, four lines 
(including the overshoot) are affect > ' . 

3. Cubic convolution can effectively remove discontinuities in the 
data due to subpixel offsets. Such offsets are introduced Into 
the Landsat MSS data by sampling delay, differences 

In scan line length, and the intersweep earth rotation effect. 

Nearest neighbor resampling cannot effectively remove thec' 
discontinuities and, as a result, photographic products gei. a ted 
from the cubic convolution data are visually much smoother than 
comparable products generated from nearest neighbor data. 

Under appropriate conditions, nearest neighbor data gives the 
appearance of containing two pixel discontinuities. 

Incorporation of a correction for these distortions is optional in 
either resampling algorithm. It was not included in the resamplings 
of data sets A, B, and C since the characteristics Investigated in 
these data sets (spreading, effect on spatial frequency, etc.) do 
not depend on subpixel effects. However, to permit a more accurate 
temporal registration, the subpixel effects were compensated for 
in data set D and E. The compensation was incorporated in both the 
cubic convolution and nearest neighbor algorithms* A comparison of 
linear ft itures In data sets D and E demonstrates the smoothing effect 
cubic convolutions can provide in photographic products. 
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4. Th« photolnterpratcr who evaluated the photographic products 
found that* compared to the original data, nearest neighbor 
TeBampllnp produces a geometric "dappling" of the finer detj>.lls 

In the data. There is not a smooth transition among neighboring 
pixels of slightly differing gray-scale values (see Figure 52) . 
Straight edges such as field patterns, runways and roads 
appear more stepped and jagged than in the original data 
(Figure 52, Figure 18). At times however, there appears to 
be an enhancement over the original data of smaller terrlan 
patterns (northeast quadrant , Figure 86) . 

The cubic convolution algorithm reduces the "dappled" appearance 
of small terrain details (Figure 63 and Figure 67) and makes 
small stream patterns appear more realistic (Figure 67) . It 
reduces the stepped appearance along the edges of fields 
(Figure 65 and Figure 69) and on airfield runways (Figure 18 
and Figure 19), and it appears to promote natural appearing 
patterns in areas with finer terrain detail (Figure 87, northeast 
quadrant) . In small, single pixel areas of high contrast rel- 
ative to their surroundings, the algorithm tends to "overshoot" 
and produce radiometric values of almost certainly misleading 
contrast (Figure 7a, 7b). 

As a result, it appears that cubic convolution is preferable 
to nearest neighbor in processing the image data to further its 
interprei-abllity . Cubic convolution does the better job of 
smoothing edges knoxm to be straight. It promotes a gradational 
gray-shade transition among smaller terrain details which Increases 
the interpreter’s ability to recognize otherwise obscure and 
subtle terrain patterns. However, it has the negative attribute 
of "overshooting" small areas of high contrast relative to the 
surroundings, and thereby promotes a misleading pattern. 

5. Cubic r .i\olution attenuates the high spatial frequencies, whereas 
nearest neighbor has essentially no effect on spatial frequencies. 
The attenuation caused by cubic convolution was nearly Identical 
in both the horizontal and vertical directions and was noticed 
equally In all cases tested. Given that the transforms which were 
compu*‘ed used 128 points, the lowest spatial frequency detected 

in the horizontal direction (excluding the 'DC' or constant term) is 

1 _ 1 

128 Sampl *.a x 57.2 Meters/Sample ” 7321.6 Meters 

an the highest spatial frequency Is 

1 _ 1 

2 Samples x 57.2 Meters/Sample * 114.4 Meters 
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within this range, the power plote ahow that cubic convolution 
begins attenuating spatial frequencies at a frequency of 
approximately 

44 1 

(128 X 57.2) Meters ■ 166.4 Meters. 

The attenuation becomes progressively more severe until the power 
In thr highest frequency (1/114.4 Meters) Is attenuated by approx- 
imately a factor of 2. 

The same occurence Is noticed In the vertical direction. In this 
case, the range of spatial frequencies Is 

1 1 1 1 
128 X 79 " 10112 Meters 2 X 79 “ 158 Meters 

The attenuation begins at a spatial frequency of approximately 

44 1 

(126 X 79) “ 229.82 Meters 

and the attenuation la again approximately a factor of tvo at 
the highest frequency. 

6. In areas of near-uniform radiance^ the differences In the output 
due to choice of resampling algorithm are minimal. Cubic con- 
volution data does have a lower standard deviation In such areas 
(and In the Image as a whole). As a result, the cubic convolution 
algorithm may be considered preferable in that It reduces the 
noise In the Input data. 

7. The multlspectral classification experiment performed as part 

of this Investigation suggests that the differences due to choice 
of resampling algorithm In the proportion estimates for various 
ground cover types are Insignificant (see Figure 70). 

Analysis of the data shows that: (a) although the differences 

in proportion eswimates are not significant, the number of pixels 
classified differently in the two resampled versions is large 
for each scene; and, (b) the differences in the two resampled 
versions are characterized by zero mean Causslan nclse. 

The interiors of the fields in cubic convolution Images are more 
homogeneous (lower standard deviations) and the transitions at 
the field boundaries are smoother as compared to the nearest 
neighbor images. Comparison of the field structure in the Images 
with the ground truth maps showed the cubic convolution images to 
be ’truer* to the scene. This is a qualitative judgment based 
on visual comparisons of field sizes and boundaries. 
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Th€ claaalflcatlon reaulta obtainad should be qualified by the 
following cons Iderat lone. 

a. The classification procedure is uaually based on a model for 
the data* The model generally used with remotely sensed 
agricultural data assumes that the data corresponding to 
each distinct ground cover type constitute a normal (mixture) 
population. This model was adopted in this study even though 
it fails the test of statistical hypothesis that the mean 
vectors of any twc fields of the same type have been drawn 
from a single normal distribution* This points to a weakness 
in the model for the data on which the classification was based 
and» therefore, a weakness in using these classification re- 
sults as a basis for evaluation of the resampling techniques. 

b* From the viewpoint of classification, the smoothing mechanism 
of cubic convolution would introduce additional spatial cor- 
relation in the data and reduce the within-field variability* 

Note that even without resampling, spatial correlation is 
introduced in the data by the sensor dynamics, atmosphere, 
and the current sampling procedures* The effect of increased 
correlation on classification is difficult to study analytically; 
however, the following qualitative remarks apply. 

o The training data, corresponding to contiguous pixels 
in a field, are assumed to be independent* Clearly, 
with cubic convolution resampling, th'^s assumption is 
harder to justify. The estimates of the class covariance 
matrices obtained as though the data were independent 
would be biased and the spread in the data would be 
underestimated * 

o The smoothing process would tend to decrease the 

^distance* (e*g., divergence, Bhattacharyya distance) 
between different classes. This would be particularly 
true for scenes with narrow fields. 

o The separability of data corresponding to different 

ground covers in a scene changes with time. The impact 
of the smoothing process of cubic convolution on class- 
ification, therefore, depends upon the time of data 
acquisition and would be greater when this separability 
is small and vice versa. 

With the training areas defined so as not to include boundary pixels, 
the nearest neighbor resampling would not affect the training sta- 
tistics and its impact on classification would be seen only in the 
boundary pixels whose identities got changed in the resampling 
process. Clearly, this impact would vary depending upon the field 
structure in the scene. 
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In sunsnary, the classification results have not provided a con* 
elusive basis for the evaluation of the two resampling techniques. 
The images generated by cubic convolution are visually of a superior 
quality and * truer* to the scene on the ground. This» however, 
does not necessarily mean that the cubic convolution images can 
be classified more accurately by a given classifier. Nor does it 
mean that a classifier can not be developed to give better class* 
ification results with the cubic convolution images as compared 
to those obtained with the nearest neighbor images > or as compared 
to the results obtained with the given classifier with the images 
generated by either resampling procedure. 

8. Regarding temporal registration, images produced by cubic 
convolution resampling have a smoother visual quality and, 
correspondingly, a lower standard deviation in near*uniform 
areas. If the images generated by a temporal registration 
correction are to be differenced and used for change detection, 
this smoothing effect of cubic convolution may be desirable, 
since it may reduce the "speckling" of the difference image. 
Alternatively, if the temporally registered Images are to be 
used as additional channels in a multlspectral classification 
experiment, results presented in the preceeding paragraph 
suggest that either resampling algorithm may work equally 
well. 
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Section 3. RECOMMENDATIONS 


1* Variations in the present cubic convolution algorithm should be 
investigated to determine whether they would lead to improved 
results. Possible variations are: 

a. Modify the cubic convolution Interpolation done in the 
horizontal direction to account for oversampling. As 
presently implemented, identical interpolations are per- 
formed in both the horizontal (oversampled) and vertical 
(non-over sampled) directions. 

b. Apply cubic convolution in the horizontal direction only. 

It is in this direction that subpixel discontinuities are 
known to exist, and, by not resampling in the vertical 
direction, the undesirable overshoot occurrences noticed 
in Data Set A may be avoided. Processing loads would be 
significantly reduced if this were done. 

2. A more complete investigation of the effects of resampling on 
multlspectral classification should be performed. The proportion 
estimates obtained with the two temporally separated images of 
this experiment (either resampled vc^rion) are significantly 
different from each other. These differences arise .ae to the 
changes with time in the spectral separability of data for dif- 
ferent ground cover types in the scene. The results obtained 
with the earlier image are far too different from the ground 
truth in each case to be of any value in a study of the Impact 

of resampling procedures. Those obtained with the second image, 
though closer to the ground truth by comparison, only constitute 
a sample of one from which to draw a conclusion. A meaningful 
conclusion may be possible from a larger study involving several 
scenes with multiple acquisitions of Imagery over each. 

3. To determine which resampling method reconstructs ”truth” more 
accurately in the photographic products, ^^truth" must be known. 
High resolution aircraft photography of an area, taken at the 
same time and in the same spectral band, may be used as a refer- 
ence. Such an evaluation should be done to determine whether 
cubic convolution is introducing patterns (e.g., moires) into the 
resampled data. 
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Figure 3. Contact print of Reston 


OR NN CC 



band 4 


band 5 


mirrors, OR,NN,CC-1800-15081-4,5 


Aberdeen Proving Grounds 


Suburban Road Network, 
Wheaton, MD 


Andrews Air Force Base 


Figure 4. 


Contact print of Aberdeen, Wheaton and Andrews AFB, 
OR, NN.CC- 1080-15192-4 






Figure 5b. lOX Enlargement of Reston Mirrors, OR- 1800- 1508 1-5 
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Figure 6a. iOX Enlargeiieiil: of Reston Mirrors, MM-1800**lS081-4 
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figure 7a. lOX Enlargefflent of Reston Mirrors, a>180<^l5081-4 



Figure 7b. lOX Enlargement of Reston Mirrors, CC-i800-lS081-5 
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Figure 12. lOX Enlargement of Aberdeen Area, NN- 1080- 15192-4 
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Figure 15. lOX Enlargement of Wheaton Area 
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lOX Enlargement of Andrews AFB, OR- 
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Figure 18. ICi.^ Enlargement of Andrews AFB, MS-1080-15192' 



Figure i9. lOX r,nlargement of Andrews AFB, CC- 1080-15192-^1 
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Figure 44. Contact print of LACIE Intensive test site, 

Hand County I2, South Dakota. Reference Image 
OR,MN,CC-2183-16433-4,5,6,? 
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Figure 70. Results of the Multispectral Classification Experiment 
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Appendix B 

Data Set A (Reston Mirrors) 
Resampled at 1:1 Pixel Aspect Ratio 



Section B.l 


SUMMARY OF EFFORT 


This appendix contains results of additional processing requested by 
personnel at the NASA Goddard Space Flight Center during their review 
of the draft of the body of this report. The additional processing 
consisted of resampling one band of the original data of the Reston 
mirrors of Data Set A to produce arrays with output pixels having a 1:1 
aspect ratio. The corresponding arrays presented earlier in this report 
have asoect ratios of 1.4:1 (approxLaate ground distance per pixel of 
57.2 meters horizontally by 79.0 meter vertically). 

The resampling to produce the desired arrays was perfcviued using both 
the nearest neighbor and cubic convolution resampling algorithms and 
was done so that output pixels represent 50.8 meters square. The results 
were plotted on film using lOOMm square spot, contact printed (Figure 
B.l), and enlarged by lOX (Figures B.2 and B.3). Pixel value print- 
outs of the area immediately surrounding the mirrors in the two 
resamplcr. versions were generated and are presented in Figures B.4 and 
B.5. 


-'O included in this appendix are results which compare the spatial 
frequency content of 1:1 data with that of 1.4:1 data. These results 
are based on power spectra of the resampled versions of the Reston 
mirrors area and were generated under funding from IBM*s Internal 
Rp«?earch and Development (IRAD) program in digital image processing. 

The method of computing the spectra was identical to that done for 
data set B in the body of this report. For all cases in this appendix, 
the spectra were averagea over 64 lines. 

Figures B.6 through B.25 present pairs of plots. The first plot of 
each pair is that of a line which passes through the lower right (south- 
eastern) mirror. The aecond plot of each pair is a power spectrum 
averaged over 64 lines of data centered on the line that is plotted. 

The plotted line should not be interpreted as being typical of all 64 
lines used to compute the spectrum. Because it does contain the mirror 
spike, the plotted line may have more power in the high frequencies. 

All plots are based on data representing the same ground area. Because 
the number of pixels representing that ground area varies with pixel 
scale and aspect ratio, the plots vary in size. One pair of plots is 
presented for horizontal lines, and one for vertical lin<s tor each 
of the following versions of the data: original, nearest neighbor 1.4:1, 

nearest neighbor 1:1, cubic convolution 1.4:1, and cubic convolution 1:1. 



Section B.2 


EVALUATION AND CONCLUSIONS 


Inspection of the generated data led to the following conclusions: 

1, The spreading and overshoot properties of the cubic convolution 
resampling algorithm noticed in the 1.4:1 data in the main repoi t 
are clearly present in the 1:1 data in this appendix. The 
spreading is seen in the southeastern mirror, which affects six 
horizontal lines in the cubic convolution 1:1 output data 
(Figure B.3), although it is present in only one horizontal 

line in the original data. The overshoot is seen in the black 
borders above and below the mirrors. The intensity values 
assigned to certain pixels surrounding the mirrors are, for 
nearest neighbor data, in the range 13-15. Due to the over- 
shoot, the same pixels in cubic convolution are in the range 1-4. 

2. The power spectrum plots show that there is significantly 
more high frequency power in the vertical direction than in 
the horizontal direction in all cases. This is true both of 
the spectra in this appendix and of those in the body of the 
report. The greater high frequency power in the vertical 
direction is to be expected because of the horizontal over- 
sampling and because of the fact that, vertically, detector-to- 
detector differences would appear as high-frequency noise. 

Because of this result, and the fact that cubic convolution 

is a mild low pass filter, it is not surprising that the atten- 
uation of the highest frequencies is greater in the ^'ertical 
direction (a factor of about 4 or 5) than it is in the horizontal 

direction (a factor of less than 2 ) when comparing the 1.4:1 

cubic convolution data to the original. 

A comparison of the spectra generated from original, nearest 
neighbor 1:1, and nearest neighbor 1.4:1 data shows that the 
power contained in a given spatial frequency is approximately 
the same in the three cases. The 1:1 data results in a better 
estimate of the spectrum. It has somewhat less power in the 
high frequencies because, with the increased sampling rate 
used to achieve the 1:1 aspect ratio, the folding frequency 
is moved out in the spectrum and the components folded back 
in have less power. 

A similar situation is noticed comparing the cubic convolution 
1.4:1 and 1:1 spectra. Lov" frequency power terms are very 
similar, but, in the 1:1 cast where the sampling was done more 
densely, the power in the high frequency terms is slightly 
lower than the power in the corresponding 1.4:1 terms due to 
the lack of additional folding. 

In a comparison of the nearest neighbor 1:1 spectra with the 
cubic convolution 1:1 spectra, the greater power in the high 
frequency terms of the former is indicative of the fact that 
nearest neighbor does not filter out the high frequency mirror 
image components of the spectrum which appear due to the 
original sampling operation. 


B.2 



Figure B.l. 
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